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The  yellow-fever  mosquito  Aedes  aegypti  is  a  major  vector  of  human  diseases,  such  as  dengue,  yellow  fever, 
chikungunya  and  West  Nile  viruses.  Chemoreceptor  organs  on  the  labella  and  tarsi  are  involved  in  hu¬ 
man  host  evaluation  and  thus  serve  as  potential  foci  for  the  disruption  of  blood  feeding  behavior.  In 
addition  to  host  detection,  these  contact  chemoreceptors  mediate  feeding,  oviposition  and  conspecific 
recognition;  however,  the  molecular  landscape  of  chemoreception  in  these  tissues  remains  mostly 
uncharacterized.  Here  we  report  the  expression  profile  of  all  putative  chemoreception  genes  in  the 
labella  and  tarsi  of  both  sexes  of  adult  Ae.  aegypti  and  discuss  their  possible  roles  in  the  physiology  and 
behavior  of  this  important  disease  vector. 
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1.  Introduction 

Aedes  aegypti  (L.)  (Diptera:  Culicidae)  mosquitoes  contribute  to 
the  spread  of  dengue,  yellow  fever,  chikungunya  and  West  Nile 
viruses  through  blood  meal-mediated  viral  transmission.  The  global 
health  risks  associated  with  this  disease  vector  are  increasing  (Bhatt 
et  al.,  2013;  Guzman  et  ah,  2010;  Barrett  and  Higgs,  2007),  as  these 
mosquitoes  thrive  in  urban  environments.  With  potential  vaccines 
still  in  development  (Wallace  et  ah,  2013),  bite  prevention  is 
essential  to  curb  the  spread  of  these  deadly  diseases. 

Disruption  of  mosquito  host-seeking  and  feeding  behavior  has 
been  achieved  through  use  of  several  repellent  chemicals  like  DEET 
and  Picaridin  (Dickens  and  Bohbot,  2013);  and  while  progress  has 
been  made  towards  understanding  the  molecular  associations  of 
these  repellents,  their  mechanisms  of  action  remains  unclear.  The 
repellency  response  appears  to  involve  multiple  chemoreception 
pathways.  Ae.  aegypti  avoidance  of  volatile  DEET  is  mediated  by  the 
Odorant  Receptor  gene  family,  a  molecular  pathway  also  involved 
in  DEET  perception  in  the  dipteran  relative  Drosophila  melanogaster 
(Ditzen  et  ah,  2008;  DeGennaro  et  ah,  2013).  Transgenic  Ae.  aegypti 
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lacking  this  olfactory  pathway  will  land  on  DEET  treated  human 
skin,  but  will  not  blood-feed  after  contact  with  DEET,  suggesting 
that  non-olfactory  molecular  pathways  of  the  tarsi,  labella  or 
elsewhere  are  involved  in  DEET  sensation  and  behavioral  avoidance 
(DeGennaro  et  ah,  2013).  Recently,  Ionotropic  Receptors  have  also 
been  shown  to  play  a  role  in  DEET  avoidance  in  D.  melanogaster 
(Kain  et  ah,  2013),  highlighting  the  need  to  investigate  the  role  of  all 
known  insect  chemoreception  gene  families  in  the  avoidance  of 
repellents  by  mosquitoes. 

Interestingly,  DEET  avoidance  by  D.  melanogaster  in  feeding  as¬ 
says  requires  at  least  three  Gustatory  Receptors  (Lee  et  ah,  2010).  In 
Ae.  aegypti,  bitter  sensing  gustatory  receptor  neurons  (GRNs)  of  the 
labella  respond  to  DEET  and  other  repellents  in  electrophysiological 
recordings  (Sanford  et  ah,  2013),  but  this  response’s  effect  on 
mosquito  behavior  is  unknown.  Thus,  putative  chemoreception 
genes  expressing  in  the  GRNs  of  the  major  gustatory  appendages 
may  serve  as  useful  targets  to  develop  novel  deterrents. 

Several  classes  of  genes  involved  in  chemoreception  have  been 
identified  in  insects:  gustatory  receptors  (Grs)  (Clyne  et  ah,  2000; 
Dunipace  et  ah,  2001;  Scott  et  ah,  2001;  Robertson  et  ah,  2003), 
odorant  receptors  (Ors)  (Clyne  et  ah,  1999;  Vosshall  et  ah,  1999; 
Robertson  et  ah,  2003),  ionotropic  receptors  (frs)  (Benton  et  ah, 
2009;  Croset  et  ah,  2010),  odorant  binding  proteins  (Obps)  (Vogt 
and  Riddiford,  1981;  McKenna  et  ah,  1994;  Galindo  and  Smith, 
2001),  sensory  neuron  membrane  proteins  ( Snmps )  (Rogers  et  ah, 
1997;  Vogt  et  ah,  1999;  Benton  et  ah,  2007;  Jin  et  ah,  2008), 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2014 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2014  to  00-00-2014 

4.  TITLE  AND  SUBTITLE 

The  genetics  of  chemoreception  in  the  labella  and  tarsi  of  Aedes  aegypti 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

United  States  Department  of  Agriculture,  Agricultural  Research 

Service, Henry  A.  Wallace  Beltsville  Agricultural  Research 

Center, Invasive  Insect  Biocontrol  and  Behavior 

Laboratory, Beltsville, MD, 20705 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

The  yellow-fever  mosquito  Aedes  aegypti  is  a  major  vector  of  human  diseases,  such  as  dengue,  yellow  fever 
chikungunya  and  West  Nile  viruses.  Chemoreceptor  organs  on  the  labella  and  tarsi  are  involved  in  human 
host  evaluation  and  thus  serve  as  potential  foci  for  the  disruption  of  blood  feeding  behavior.  In  addition  to 
host  detection,  these  contact  chemoreceptors  mediate  feeding,  oviposition  and  conspecific  recognition; 
however,  the  molecular  landscape  of  chemoreception  in  these  tissues  remains  mostly  uncharacterized.  Here 
we  report  the  expression  profile  of  all  putative  chemoreception  genes  in  the  labella  and  tarsi  of  both  sexes 
of  adult  Ae.  aegypti  and  discuss  their  possible  roles  in  the  physiology  and  behavior  of  this  important 
disease  vector. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

9 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


J.T.  Sparks  et  at.  /  Insect  Biochemistry  and  Molecular  Biology  48  (2014)  8-16 


9 


pickpocket  channels  ( ppks )  or  DEG/ENaCs  (Liu  et  al„  2003),  tran¬ 
sient  receptor  potential  channels  {Trps)  (Al-Anzi  et  al„  2006;  Kwon 
et  al„  2010;  Kim  et  al.,  2010;  Kang  et  al.,  2010)  and  CheA/Bs  (Xu 
et  al.,  2002;  Park  et  al.,  2006).  These  genes  are  typically  expressed 
in  chemosensory  neurons  or  support  cells  associated  with  pe¬ 
ripheral  hair-like  sensory  organs  called  sensilla.  These  chemo¬ 
sensory  neurons  discriminate  between  diverse  chemicals,  often 
informing  mosquito  behaviors  such  as  host-seeking,  feeding, 
oviposition  and  mating. 

There  are  few  comprehensive  reports  of  chemoreception  gene 
expression  in  specific  appendages  in  dipteran  species;  confirmation 
of  expression  in  D.  melanogaster  has  been  shown  through  non- 
quantitative  visual  representations  using  promoter-driven  re¬ 
porters  or  amplified  PCR  products,  with  these  reports  being  orga¬ 
nized  by  chemoreception  gene  family  rather  than  global  expression 
profile  (Vosshall  et  al.,  1999;  Galindo  and  Smith,  2001 ;  Couto  et  al., 
2005;  Benton  et  al„  2009;  Isono  and  Morita,  2010).  In  Ae.  aegypti, 
characterization  of  Or  expression  has  been  conducted  non- 
quantitatively  by  PCR  in  all  head  appendages  (Bohbot  et  al., 
2007),  and  three  Grs  have  been  studied  in  the  maxillary  palps 
(Erdelyan  et  al.,  2012;  Bohbot  et  al.,  2013).  We  previously  reported 
the  expression  profile  of  Grs  in  the  labella  and  tarsi  of  male  and 
female  Ae.  aegypti  (Sparks  et  al.,  2013).  Here  we  extend  our  survey 
to  the  expression  of  other  putative  chemoreception  genes  by 
RNA-seq  in  these  tissues,  thus  expanding  our  knowledge  of  the 
appendage-specific  molecular  components  potentially  involved  in 
chemical  attraction  and  avoidance.  Expressed  putative  chemore¬ 
ception  genes  in  male  and  female  mosquitoes  may  be  targeted  in 
future  experiments  aimed  at  altering  their  behavior  and  ultimately 
disrupting  transmission  of  harmful  viruses.  Furthermore,  we 
assessed  the  significance  of  this  expression  for  a  small  set  of 
chemoreception  genes  by  comparing  RT-qPCR  estimates  of  ex¬ 
pression  between  female  labella  and  female  carcass  tissue  samples. 

2.  Materials  and  methods 

2.3.  Animal  rearing 

Ae.  aegypti  eggs  (Orlando  strain)  were  obtained  from  the  Center 
for  Medical  and  Veterinary  Entomology,  USDA,  ARS  in  Gainesville, 
FL,  USA.  Larvae  were  reared  at  25  °C  (12-hL:12-hD)  and  fed  with 
ground  TetraMin®  fish  food.  Unsexed  pupae  were  hand-collected 
daily  and  transferred  to  plastic  dishes  (9  cm  x  5.5  cm)  inside 
small  containment  buckets,  thus  establishing  24-h  age  groups. 
Greater  than  95%  of  adults  emerged  2  days  post-pupation,  after 
which  all  remaining  pupae  were  removed  from  containment 
buckets.  Adult  mosquitoes  were  fed  with  a  10%  sucrose  solution 
and  maintained  in  an  environmental  chamber  at  27  °C  and  70% 
relative  humidity  under  the  same  photoperiod  as  larvae.  Tissues 
used  in  our  studies  were  collected  during  the  photophase  from 
adult  mosquitoes  6—7  days  old. 

2.2.  RNA  isolation  and  sequencing 

For  RNA  sequencing,  paired  labella  from  500  males  or  500  fe¬ 
males  were  carefully  dissected  to  limit  inclusion  of  other  adjacent 
proboscis  tissues.  Samples  from  legs  were  comprised  of  pro-,  meso- 
,  or  metathoracic  tarsal  segments  of  400  males  or  400  females. 
Dissected  tissues  were  immediately  stored  on  dry  ice  and  me¬ 
chanically  ground  in  TRIzol®  (Life  Technologies,  Carlsbad,  CA,  USA). 
Total  RNA  was  isolated  by  RNeasy®  Plus  Mini  Kit  (Qiagen,  Valencia, 
CA,  USA),  quantified  on  a  Nanodrop  ND-1000  spectrophotometer 
(Nano  Drop  Products,  Wilmington,  DE,  USA),  and  sent  to  the  Ge¬ 
nomics  Services  Lab  at  the  Hudson  Alpha  Institute  for  Biotech¬ 
nology  (Huntsville,  AL).  Messenger  RNA  isolation  and  subsequent 


cDNA  synthesis  were  completed  using  NEBNext®  reagents  (NEB, 
Ipswich,  MA,  USA)  and  standard  protocols  with  custom  GSL 
adaptors.  Complementary  DNA  libraries  corresponding  to  distinct 
tissues  were  sequenced  on  an  Illumina  HiSeq2000  to  generate  25 
million  50  base  pair,  paired-end  reads  per  sample. 

2.3.  Analysis  of  annotated  and  unannotated  chemosensory  genes 

Reference  genome  and  annotations  for  Ae.  aegypti  (AaegL1.3) 
were  downloaded  from  VectorBase  (http://aaegypti.vectorbase. 
org/GetData/Downloads/).  Output  Fastq  Illumina  files  were  map¬ 
ped  to  the  reference  genome  with  TopHat  (Trapnell  et  al.,  2009). 
The  unambiguous  sequence  alignment  files  were  uploaded  into  the 
Avadis  NGS  software  (Strand  Scientific  Intelligence,  CA,  USA), 
where  quantification  and  normalization  were  performed.  Prior  to 
quantification  using  the  ‘Deseq’  normalization  method,  the  read  list 
was  filtered  to  remove  duplicate,  single-end,  mate-filtered,  mate¬ 
missing,  one-mate  flip,  both-mate  flip,  and  unaligned  reads.  Read 
quality  metric  values  were  as  follows:  Quality  threshold  >30,  N’s 
allowed  in  read  <0,  Alignment  score  threshold  >95,  Mapping 
quality  threshold  >40.  Transcript  expression  levels  for  all  genes 
are  reported  in  values  of  Reads  Per  Kilobase  per  Million  reads 
mapped  (RPKM).  RPKM  values  represent  a  quantitative  measure  of 
the  number  of  corresponding  50  bp  sequence  reads  (sequenced  in 
both  directions)  for  a  given  gene.  We  assigned  no  specific  RPKM 
threshold  for  functional  expression  vs.  background  “noise.” 

2.4.  Quantitative  RT-PCR 

Thirteen  chemosensory  genes  and  one  housekeeping  gene  were 
selected  for  qPCR  analysis  to  evaluate  gene  expression  over  a  dy¬ 
namic  range,  both  in  copy  number  and  presumed  chemosensory 
gene  function.  Primer  pairs  were  designed  for  each  target  gene  to 
amplify  a  specific  100—180  base  pair  PCR  product  (Primer-BLAST 
Primer  Designing  tool,  NCBI).  At  least  one  primer  per  set  spans  an 
exon  boundary  to  exclude  non-specific  gDNA  amplification. 

Statistically  supported  RT-qPCR  validation  of  RNA  sequencing  of 
labella  and  all  tarsal  types  of  both  sexes  was  previously  reported  in 
Sparks  et  al.  (2013).  Here  we  dissected  150  paired  female  labella 
and  10  female  carcasses  (thorax,  halteres  and  abdomen)  for  RNA 
extraction,  repeating  each  tissue  once.  Total  RNA  was  isolated  from 
all  frozen  tissue  samples  as  previously  described.  cDNA  was  syn¬ 
thesized  using  Superscript®  III  First-Strand  Synthesis  Supermix  for 
qRT-PCR  (Life  Technologies,  Carlsbad,  CA,  USA).  PCR  products  were 
directly  sequenced  to  confirm  amplicon  identity  (data  not  shown) 
(Macrogen,  Rockville,  MD,  USA).  RT-qPCR  was  subsequently  per¬ 
formed  on  each  target  gene  using  KiCqStart®  SYBR®  Green  qPCR 
ReadyMix™  iQ  (Sigma— Aldrich,  St.  Louis,  MO,  USA)  and  an  iCycler 
iQ™  Real-Time  PCR  Detection  System  (Bio-Rad,  Hercules,  California, 
USA).  All  Ct  values  were  calculated  by  Bio-Rad  iQ5  Optical  System 
Software  (Bio-Rad,  Hercules,  California,  USA).  Reactions  were  per¬ 
formed  in  technical  triplicate  20  pL  volumes.  Three-step  cycles  plus 
melt  curves  were  used  for  each  reaction,  using  an  annealing  tem¬ 
perature  of  56  °C  for  all  primers.  Efficiencies  for  each  primer  set 
were  calculated  from  the  slope  of  the  standard  curve  using  the 
formula  E  =  10'  1  />|opeJ  ( pfaffl ,  2001;  Rasmussen,  2001).  Primer 
efficiencies  are  based  on  three  1:10  serial  dilutions  of  cDNA  tem¬ 
plate  used  in  side-by-side  technical  triplicate  reactions.  Efficiencies 
are  listed  in  Table  SI. 

Relative  gene  quantification  was  calculated  as  E targelt a rge c  1 — c c 
[reference])  for  eacj1  target  gene  (14  total)  and  averaged  for  each 
replicate  (6  total  for  each  tissue  sample,  representing  two  biological 
replicates).  Ae.  aegypti  housekeeping  gene  Lysosomal  Aspartic 
Protease  (Vectorbase  ID:  AAEL006169)  was  used  to  normalize  Ct 
values  between  tissue  types.  ‘Times-enrichment’  was  calculated  as 
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Fig.  1.  Prominent  chemoreception  genes  of  the  labella.  Each  colored  cell  represents  the 
average  RPKM  value  between  male  and  female  labellar  samples  (tissue  collected  in 
blue).  Only  annotated  chemoreception  genes  yielding  averages  above  3  RPKM  are 
listed  except  for  Obps,  which  are  at  least  100  RPKM.  Individual  heat  map  scales  are 
shown  to  highlight  the  most  abundantly  expressing  genes  within  each  gene  family. 
(For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred 
to  the  web  version  of  this  article.) 


the  ratio  of  female  labella  transcript  abundance  over  female  carcass 
transcript  abundance. 

3.  Results 

Here  we  report  the  labellar  and  tarsal  expression  of  genes 
associated  with  chemoreception  that  potentially  affects  behaviors 
such  as  host-seeking,  feeding,  oviposition  and  mating.  Gene 
expression  is  reported  as  RPKM,  which  denotes  relative  transcript 
abundance  for  each  putative  chemoreception  gene.  For  a  small  se¬ 
lection  of  chemoreception  genes,  we  also  quantify  relative  abun¬ 
dance  in  the  labella  with  respect  to  carcass.  We  highlight  possible 
functions  for  these  genes  drawing  on  functionally  characterized 
insect  homologs. 
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Fig.  2.  Prominent  chemoreception  genes  of  the  pro-  (A),  meso-  (B)  and  meta-thoracic 
(C)  tarsi.  Each  colored  cell  represents  the  average  RPKM  value  between  male  and  fe¬ 
male  tarsal  samples  (tissue  collected  in  blue).  Only  annotated  chemoreception  genes 
yielding  averages  above  3  RPKM  are  listed  except  for  Obps,  which  are  at  least  1000 
RPKM.  A  full  list  of  annotated  chemoreception  gene  expression  is  available  in  the 
supplementary  information.  Individual  heat  map  scales  are  shown  to  highlight  the 
most  abundantly  expressing  genes  within  each  gene  family  for  the  three  tissue  types. 
(For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred 
to  the  web  version  of  this  article.) 


3.1.  Odorant  Receptors  (Ors) 

Insect  ORs  function  as  heterodimers  to  form  ligand  gated  ion 
channels  (Sato  et  al„  2008;  Wicher  et  al.,  2008),  with  single  ORs 
conferring  odor  sensitivity  to  individual  odorant  receptor  neurons 
(ORNs).  Our  survey  showed  expression  of  the  ubiquitous  co¬ 
receptor  AaegOrco  and  several  putative  ligand-binding  AaegOrs  in 
male  and  female  labella  (Fig.  1  and  SI);  RPKM  values  for  AaegOrco 
were  similar  for  both  sexes  (123.8,  male  labella;  112.2,  female 
labella).  Six  other  AaegOrs  showed  RPKM  values  above  10.0  ( Aae - 
gOr4,  AaegOr6,  AaegOrl5,  AaegOrl8,  AaegOr46  and  AaegOr47 )  and 
three  Ors  showed  intermediate  RPKM  values  between  3.3  and  8.4 
(AaegOr2,  AaegOr81  and  AaegOrl07)  supporting  previous  reports  of 
Or  expression  in  the  proboscis  of  mosquitoes  (Pitts  et  al.,  2004; 
Melo  et  al.,  2004;  Kwon  et  al.,  2006;  Bohbot  et  al.,  2007).  Inter¬ 
estingly,  there  was  low-level  (RPKM  between  2.2  and  5.5)  expres¬ 
sion  of  a  few  putative  ligand-binding  Ors  in  the  tarsi  ( AaegOr2 , 
AaegOr52,  AaegOr62  and  AaegOrl07),  but  no  Oreo  expression  above 
0.5  RPKM  (Fig.  2  and  SI).  As  insect  ORs  likely  ubiquitously  require 
ORCO  co-expression  for  functional  activation  and  dendritic  traf¬ 
ficking  (Larsson  et  al.,  2004;  Jones  et  al.,  2005;  Benton  et  al.,  2006), 
it  is  unlikely  AaegOr  expression  in  the  tarsi  significantly  contributes 
to  chemoreception  in  this  context. 

Two  AaegOrs  that  express  in  labella  of  Ae.  aegypti  are  highly 
similar  in  sequence  to  an  Or  from  the  malaria  mosquito  Anopheles 


gambiae.  AaegOr4  and  AaegOrG  are  closely  related  to  the  AgamOr6 
(Bohbot  et  al.,  2007),  which  is  tuned  to  detect  several  naturally 
occurring  aromatics  like  acetophenone  and  heterocyclic  com¬ 
pounds  (Carey  et  al.,  2010)  and  is  expressed  in  the  proboscis  of 
Anopheles  gambiae  (Kwon  et  al.,  2006).  Thus,  AaegOR4  and  AaegOR6 
may  be  sensitive  to  similar  classes  of  volatile  chemicals  in  the 
context  of  the  labellum.  The  significance  of  these  chemicals  for 
mosquito  ecology  is  unclear. 

While  the  closely  related  (49%  identical,  69%  similar)  AaegOr46 
and  AaegOr47  receptors  share  no  homologs  with  greater  than  25% 
identity  in  A.  gambiae  (Bohbot  et  al.,  2007),  of  the  four  most  closely 
related  AgamOrs  to  this  Ae.  aegypti  Or  subfamily  (100  bootstrap 
support,  1000  iterations),  two  (Again OR56  and  AgamOR57)  also 
show  expression  in  A.  gambiae  whole  proboscis  (Kwon  et  al.,  2006). 
Whether  the  conservation  of  proboscis  expression  between  these 
mosquito  ORs  is  related  to  function  requires  further  investigation. 

3.2.  Ionotropic  Receptors  (Irs) 

IRs  are  ligand-binding  chemoreceptors  in  insects  that  have 
likely  evolved  from  animal  ionotropic  glutamate  receptors  (Croset 
et  al.,  2010),  which  are  a  conserved  family  of  ligand-gated  ion 
channels  associated  synaptic  communication  in  both  eukaryotes 
and  prokaryotes.  Since  their  identification  (Benton  et  al.,  2009), 
IRs  have  been  primarily  described  as  olfactory  receptors  in 
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D.  melanogaster  (Abuin  et  al.,  2011 ;  Ai  et  al.,  2013),  although  other 
functional  roles  have  been  assigned  to  IRs  expressing  in  non- 
antennal  tissues  (Senthilan  et  al.,  2012;  Zhang  et  al.,  2013). 

Thirty-six  putative  Irs  show  expression  levels  above  3.0  RPKM  in 
the  labella  and  tarsi  of  Ae.  aegypti  (Figs.  1  and  2  and  SI),  which 
strongly  suggests  their  involvement  in  chemoreception  in  these 
appendages.  Aaeglr25a,  which  has  clear  orthologs  in  all  pro- 
tostomes  (Croset  et  al.,  2010)  and  is  the  most  similar  1R  to  iGluRs, 
showed  the  highest  RPKM  values  for  labella,  pro-  and  meso- 
thoracic  tarsi,  with  RPKM  values  over  300  for  male  and  female 
labella  (Figs.  1  and  2).  Both  DmelIR25a  expression  (Benton  et  al., 
2009)  and  function  (Abuin  et  al.,  2011)  in  ORNs  innervating  many 
coeloconic  sensilla  on  the  antenna  of  D.  melanogaster  strongly 
suggest  lR25a  is  a  co-receptor  for  other  ligand-binding  IRs  (Benton 
et  al.,  2009).  This  function  is  consistent  with  our  expression  data,  as 
Aaeglr25a  is  the  most  abundantly  expressed  Ir  in  our  survey. 

Other  conserved  putative  Aaeglrs  (Rytz  et  al.,  2013)  showing 
expression  in  labella  and/or  tarsi  are:  Aaeglr8a,  Aaeglr76b,  the 
Aaeglr4l  subfamily,  and  two  members  of  the  Aaeglr75  subfamily 
(Figs.  1  and  2).  Aaeglr8a  shows  RPKM  values  between  5.9  and  9.6  in 
male  and  female  tarsi  (Fig.  2),  much  lower  than  Aaeglr25a.  Like 
AaegIR25a,  AaeglR8a  may  function  as  an  IR  co-receptor  based  on 
data  from  its  ortholog  DmelIR8a  (Benton  et  al.,  2009;  Abuin  et  al., 
2011).  Aaeglr76b  is  abundantly  expressed  in  the  labella  (RPKM 
reaching  153.3  in  males)  and  prothoracic  tarsi  of  males  and  females 
and  less  so  in  meso-  and  meta-thoracic  tarsi.  Aaeglr76b  is  orthol- 
ogous  to  Dmellr76b,  which  is  involved  in  salt  reception  in 
D.  melanogaster  GRNs  (Zhang  et  al.,  2013),  but  is  also  expressed  in 
all  classes  of  coeloconic  sensilla  on  the  antenna  (Benton  et  al., 
2009).  The  Aaeglr4l  subfamily  shows  RPKM  values  as  high  as  42.4 
for  Aaeglr4le,  and  most  of  these  thirteen  Irs  show  expression  above 
10.0  RPKM  in  male  and  female  labella.  Though  the  Dmellr4la/Aae- 
glr41  group  appears  to  be  conserved  throughout  several  insect  or¬ 
ders  (Rytz  et  al.,  2013),  this  is  only  the  second  report  of  its 
expression  (Leal  et  al.,  2013).  The  Ir75  subfamily  in  Ae.  aegypti  has 
expanded  to  ten  members;  two  of  these  Irs,  Aaeglr75k.2  and  Aae- 
glr75l,  are  expressed  in  labella  and  tarsi  of  males  and  females. 
Dmel IR75  is  likely  involved  in  propionic  acid  reception  in  the  an¬ 
tenna  of  D.  melanogaster  (Yao  et  al.,  2005;  Abuin  et  al.,  2011).  Pro¬ 
pionic  acid  and  other  carboxylic  acids  are  mosquito  attractants  in 
the  context  of  odor  blends  containing  ammonia  and  lactic  acid 
(Smallegange  et  al.,  2009).  Thus,  AaegIR75k.2  and  AaeglR751  may  be 
involved  in  the  detection  of  attractants  emitted  from  human  skin, 
which  include  several  aliphatic  carboxylic  acids  (Bosch  et  al.,  2000). 

In  addition  to  these  conserved  Irs,  eleven  of  the  nineteen  Aaeglrs 
(with  accession  numbers)  associated  with  the  dipteran  Ir7  sub¬ 
family  (Croset  et  al.,  2010)  are  expressed  at  RPKM  levels  above 
4.0  in  male  and  female  labella  (Fig.  1 ).  This  expression  is  consistent 
with  the  observation  that  Ir7  subfamily  member  Dmellr7a  is 
expressed  in  both  adult  and  larval  gustatory  neurons  in 
D.  melanogaster  (Croset  et  al.,  2010).  Thus,  it  likely  that  at  least  some 
of  the  expressed  AaeglRs  are  involved  in  chemical  reception  in  the 
labellum. 

Other  putative  Irs  showing  RPKM  values  higher  than  3.0  in 
labella  and/or  tarsi  are  AaegIrI20,  AaegIrlII.2,  Aaeglrl03,  Aae- 
glrl04.2,  Aaeglrl04.3,  and  five  unnamed  Irs  (Figs.  1  and  2).  Among 
dipteran  genomes,  Aaeglrl20  is  most  closely  related  to  Cquilrl20 
(Croset  et  al.,  2010),  which  is  more  abundant  in  the  metathoracic 
tarsi  than  the  antennae  of  Culex  quinquefasciatus  (Leal  et  al.,  2013). 
This  is  consistent  with  tarsal  expression  levels  for  Aaeglrl20 
observed  here  (Fig.  2). 

Genes  specific  to  or  enriched  in  chemoreception  organs  of  fe¬ 
male  Ae.  aegypti  are  candidates  for  involvement  in  blood  host 
location  and  recognition.  There  were  no  obvious  sex  biases  among 
AaegOrs  expressed  in  labella,  however  a  few  Aaeglrs  demonstrated 


sex  biases.  AaegIrI03,  AaegIrI04.2  and  Aaeglrl04.3  are  most  abun¬ 
dant  (>4  times  enriched)  in  female  pro-  and  meso-thoracic  tarsi 
and  are  not  expressed  in  labella  of  either  sex  (Figure  SI).  The  tarsi 
represent  the  initial  point  of  contact  between  feeding  female 
mosquitoes  and  their  hosts  (Clements,  1992);  thus,  tarsi  may  house 
female-specific  chemoreceptors  that  first  assess  the  non-volatile 
cues  associated  with  potential  blood  hosts.  Alternatively,  these 
female  tarsal  enriched  IRs  may  be  involved  in  the  recognition 
of  conspecific  males,  a  function  of  the  tarsi  of  D.  melanogaster 
(Greenspan  and  Ferveur,  2000). 

3.3.  Sensory  neuron  membrane  proteins  (Snmps) 

Snmps  are  members  of  a  gene  family  characterized  by  human 
CD36,  having  diverse  functions  including  transport  of  fatty  acids 
and  interactions  with  other  proteins  in  celkcell  recognition  (Vogt, 
2003;  Silverstein  and  Febbraio,  2009).  SNMPs  are  membrane  pro¬ 
teins  observed  to  associate  with  chemosensory  sensilla  in  insects, 
and  are  conserved  throughout  four  orders  of  the  holometabolous 
lineage  (Vogt  et  al.,  2009).  In  D.  melanogaster,  SNMP1  is  essential  for 
the  detection  of  the  pheromone  cis-vaccenyl  acetate  (Benton  et  al., 
2007;  Jin  et  al.,  2008).  Here  we  focus  on  the  two  putative  Snmps 
identified  in  Ae.  aegypti,  Snmpl  and  Snmp2  (Vogt  et  al.,  2009). 

AaegSnmp 2  is  abundantly  expressed  in  the  labella  and  tarsi  of 
both  sexes,  although  RPKM  values  are  somewhat  higher  in  males 
(RPKM  range  140.8—168.9  for  male  tissue,  81.9—104.3  for  female 
tissue;  Figs.  1  and  2).  This  is  consistent  with  reported  Snmp2 
expression  in  leg  and  gustatory  tissue  from  the  dipteran 
D.  melanogaster  (Benton  et  al.,  2007;  Vogt  et  al.,  2009)  as  well  as  the 
lepidopterans  Cnaphalocrocis  medinalis  (Liu  et  al.,  2013)  and  Agrotis 
ipsilon  (Gu  et  al.,  2013).  As  yet,  no  function  has  been  attributed  to 
SNMP2  in  insects. 

Conversely,  AaegSnmp  I  showed  almost  no  expression  in  labella 
or  tarsi  (maximum  RPKM  of  0.1  in  female  labella  and  metathoracic 
tarsi;  Figure  SI),  an  unexpected  result  considering  a  previous  non- 
quantitative  report  of  leg  expression  in  Ae.  aegypti  (Vogt  et  al., 
2009).  It  is  possible  that  AaegSnmpl  transcript  was  represented  in 
femoral  or  tibial  tissues  in  this  instance,  as  these  leg  segments  were 
not  included  in  our  survey. 

Snmpl  homologs  are  expressed  more  broadly  in  other  insect 
species.  Relatively  low  expression  in  bodies  and  legs  of 
D.  melanogaster  has  been  reported  (Benton  et  al.,  2007;  Vogt  et  al., 
2009).  Snmpl  expression  in  lepidopteran  species  varies  from 
widespread,  including  legs,  in  Plutella  xylostella  (Li  and  Qin,  2011), 
C.  medinalis  (Liu  et  al.,  2013)  and  A.  ipsilon  (Gu  et  al.,  2013),  to 
antennal-specific  in  other  species  including  A.  gambiae  (Benton 
et  al,  2007),  Manduca  sexta  (Rogers  et  al.,  2001 )  and  Amyelois  tran- 
sitella  (Leal  et  al.,  2009). 

3.4.  Ion  channels  associated  with  chemoreception  (Trps  and  Ppks) 

TRP  channels  belong  to  a  highly  conserved  family  of  ion  chan¬ 
nels  that  express  and  function  in  multiple  types  of  sensory  neurons 
in  insects  (Fowler  and  Montell,  2013).  Here  we  focus  on  potential 
roles  in  chemoreception,  although  these  membrane  proteins  are 
also  involved  in  photo-,  thermo-,  and  mechanoreception  in  insects. 
Five  putative  AaegTrp  channels  (Bohbot  et  al.,  2014)  show  expres¬ 
sion  above  3.6  RPKM  in  male  and  female  labella  ( Trpal ,  Painless, 
Trpml,  one  of  the  two  Nompc  homolgues,  and  Trpy)  (Fig.  1).  All  but 
Trpy  also  express  above  4.7  RPKM  in  at  least  one  tarsal  sample 
(Fig.  2).  Lower  RPKM  values  may  indicate  functional  significance  for 
other  putative  AaegTrp  channels  (Figure  SI ). 

AaegTRPAl,  AaegPainless  and  AaegTRPy  are  the  most  likely 
TRPs  to  contribute  to  chemoreception  in  labella  and  tarsi  as  ho¬ 
mologs  of  these  channels  affect  chemoreception  in  other  insects.  In 
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D.  melanogaster,  neuronal  responses  to  citronellal  in  the  antenna 
are  mediated  by  TRPA1  (Kwon  et  al.,  2010).  While  citronellal  did  not 
directly  activate  D.  melanogaster  TRPA1  in  the  aforementioned 
study,  it  did  directly  activate  mosquito  A.  gambiae  TRPA1.  Thus, 
TRPA1  may  contribute  in  some  way  to  the  reception  of  repellents 
by  GRNs  in  the  context  of  the  labella  and  tarsi  of  Ae.  aegypti. 
Similarly,  noxious  chemical  electrophiles  have  been  shown  to  deter 
feeding  responses  in  D.  melanogaster  and  evoke  Dme/TRPAl-  and 
AgamTRPAl-mediated  responses  when  they  are  expressed  in  a 
heterologous  system  (Kang  et  al.,  2010).  Furthermore,  Dme/TRPAl 
appears  to  be  involved  in  gustatory  reception  of  aristolochic  acid 
but  not  broadly  for  other  bitter  compounds  (Kim  et  al.,  2010).  In  the 
lepidopteran  Manduca  sexta,  TRPA1  mediates  gustatory  responses 
to  aristolochic  acid  and  integrates  information  about  temperature 
(Afroz  et  al.,  2013),  possibly  a  result  of  altered  thermosensitivity 
through  the  evolution  of  novel  isoforms  to  gain  or  maintain 
chemical  sensitivity  (Kang  et  al.,  2012). 

Dme/Painless  is  required  for  the  gustatory-mediated  avoidance 
of  isothiocyanate  (Al-Anzi  et  al.,  2006);  thus  AaegPainless  (the 
closest  Ae.  aegypti  homolog,  Bohbot  et  al.,  2014)  expression  in  the 
labella  may  indicate  a  similar  sensitivity  in  mosquitoes. 

Dme/TRPL,  a  member  of  the  TRPC  subfamily  associated  with 
phototransduction  (Fowler  and  Montell,  2013),  is  also  involved  in 
the  feeding  adaptation  to  camphor  in  D.  melanogaster  (Zhang  et  al., 
2013b).  Dme/TRPL  is  closely  related  to  AaegTRPy  (Bohbot  et  al., 
2014);  therefore,  AaegTrpy  expression  in  the  labella  may  indicate 
involvement  in  gustatory  perception. 

Ppk  channels  are  members  of  the  diverse  animal-specific  DEG/ 
ENaC  family  of  ion  channels.  While  poorly  understood  in  insects  in 
general,  D.  melanogaster  Ppks  are  spatially  and  temporally  diverse 
in  their  expression  pattern  (Zelle  et  al.,  2013),  and  are  involved  in 
the  gustatory  detection  of  salt  (Liu  et  al.,  2003)  and  mating  cues  (Lin 
et  al.,  2005;  Thistle  et  al.,  2012;  Lu  et  al.,  2012;  Toda  et  al.,  2012; 
Starostina  et  al.,  2012).  We  observed  moderate  expression  (RPKM 
range  3—27)  of  10  Aaegppks  in  labella  and  7  Aaegppks  (RPKM  range 
3—11)  in  tarsi  (Figure  SI).  Of  these  genes,  only  2  Aaegppks  are 
expressed  in  both  tissue  types.  We  also  detected  low  Aaegppk 
expression  (RPKM  1—3)  in  multiple  instances  (Figure  SI),  which 
may  also  indicate  functional  significance  in  a  restricted  number  of 
cells. 

3.5.  Small,  soluble  proteins  associated  with  chemoreception  (Obps 
and  CheA/Bs) 

Insect  OBPs  are  small  soluble  proteins,  many  of  which  are 
excreted  into  the  lymph  surrounding  chemoreceptor  neurons  and 
influence  the  detection  of  both  volatile  (Vogt  and  Riddiford,  1981; 
Vogt  et  al.,  1999)  and  non-volatile  compounds  (Jeong  et  al.,  2013; 
Swarup  et  al.,  2013).  The  Ae.  aegypti  genome  contains  111  Obp 
genes  subdivided  into  three  subfamilies  that  have  been  defined 
based  on  bioinformatics  data  (Zhou  et  al.,  2008;  Manoharan  et  al., 
2013).  Classic  OBPs  have  a  conserved  motif  of  6  cysteine  residues. 
Plus-C  OBPs  possess  additional  cysteine  residues,  and  two-domain 
(also  referred  to  as  “atypical")  OBPs  are  the  largest  OBP  subfamily. 

AaegObps  showed  the  highest  level  of  transcript  abundance  in 
our  survey  (RPKM  range  0—45,000;  Fig.  3).  While  there  are  some 
differences  in  expression  levels  between  sexes  (Figure  SI),  sex- 
biased  expression  of  Obps  likely  does  not  account  for  physiolog¬ 
ical  or  behavioral  differences  between  male  and  female  Ae.  aegypti. 
The  most  abundantly  expressed  AaegObps  in  the  labella  of  both 
sexes  are  Obp57,  Obp69,  Obp42,  Obp63,  Obp24,  Obp25,  Obp72  and 
Obpll  (RPKM  range  3500—45,000;  Figs.  1  and  3).  Some  Obps,  like 
Obp 5,  Obpll,  Obp35  and  Obp81,  were  more  abundant  in  specific 
tarsal  types  than  in  labella,  often  a  feature  shared  between  sexes  for 
specific  genes.  In  all,  18  Obps  showed  RPKM  values  greater  than  100 


in  labella  and  16  showed  RPKM  values  above  100  in  at  least  one 
tarsal  tissue  (Figure  SI).  The  two-domain  Obps  are  generally  not 
expressed  in  labella  or  tarsi  (Fig.  3),  possibly  highlighting  a  func¬ 
tional  difference  in  this  sub-family  of  Obps. 

Functional  data  on  mosquito  OBPs  is  lacking.  The  AaegObpl 
homolog  An.  gambiae  Obpl  has  been  reported  to  specifically  bind 
to  the  oviposition  attractant  indole  (Biessmann  et  al.,  2010).  Aae¬ 
gObpl  homolog  C.  quinquefasciatus  Obpl  has  been  implicated  in  the 
sensitivity  of  female  mosquitoes  to  several  oviposition  cues 
(Pelletier  et  al.,  2010).  We  observed  relatively  little  AaegObpl 
expression  in  the  tissues  we  surveyed  (Fig.  3).  The  knockdown  of 
two  Ae.  albopictus  Obps  ( AalbObp37  and  AalbObp39 )  affected 
antennal  responses  to  indole  (Deng  et  al.,  2013).  AaegObp39,  the 
closest  homolog  of  AalbObp39,  showed  moderate  expression  in  all 
tissues  we  surveyed  (Fig.  3). 

Insect  CheAs  and  CheBs  represent  two  separate  gene  families 
that  encode  small,  membrane-bound  proteins  that  are  excreted  in 
some  cases  (Xu  et  al.,  2002;  Ben-Shahar  et  al.,  2010).  These  genes  are 
distinct  from  Obps  though  they  too  are  usually  associated  with 
chemoreceptor  organs.  Some  CheA  and  CheB  proteins  are  expressed 
in  sex-specific  patterns  and  are  thought  to  mediate  pheromone 
reception  in  sex-specific  tarsal  sensilla  of  D.  melanogaster  (Xu  et  al., 
2002;  Park  et  al.,  2006).  The  molecular  mechanism  of  CheA  and 
CheB  function  remains  unknown,  but  may  involve  interaction  with 
other  membrane  bound  receptors  or  Ppk  channels  (Ben-Shahar 
et  al.,  2010). 

As  in  D.  melanogaster,  most  AaegCheA  and  AaegCheB  genes  ex¬ 
press  in  tarsi  in  a  sex-specific  manner.  These  genes  likely  influence 
mating  behavior  by  selectively  binding  to  and  transporting  non¬ 
volatile  pheromones  required  to  initiate  successful  copulation  (Xu 
et  al.,  2002;  Park  et  al.,  2006).  Of  ten  putative  AaegCheA  genes 
(gene  curation,  Bohbot  et  al.,  2014),  four  are  highly  enriched  in  male 
tarsi  with  respect  to  female  tarsi  and  one  is  highly  enriched  in 
female  tarsi  with  respect  to  male  tarsi  (Figure  SI).  The  remaining 
five  AaegCheA  genes  are  abundant  in  tarsi  of  both  sexes.  All  three 
putative  AaegCheB  genes  are  enriched  in  male  tarsi  with  respect  to 
female  tarsi,  or  the  labella  of  both  sexes  (Figure  SI). 

3.6.  Enrichment  of  chemoreception  genes  in  female  labella 

Enrichment  of  gene  transcripts  in  the  labellum  with  respect  to 
carcass  supports  their  involvement  in  chemoreception.  We 
compared  the  expression  of  14  genes  through  RT-qPCR  of  female 
labella-  and  female  carcass-derived  cDNA  libraries  (Fig.  4);  relative 
abundance  of  each  gene  is  presented  as  the  ratio  of  labella  to 
carcass  expression.  As  expected,  the  two  genes  showing  near  zero 
RPKM  values  for  female  labella  RNA-seq  samples  (AaegGRl,  0.2; 
AaegGR76,  0.0)  also  showed  the  lowest  enrichment  in  female 
labella  tissue  by  RT-qPCR  (AaegGRl,  1.1-times;  AaegGR76,  4.8- 
times).  All  other  putative  chemoreception  genes  in  this  target 
group  showed  relatively  high  enrichment,  between  13  and  15,000- 
times,  in  female  labella  and  correspond  to  relatively  high  RPKM 
values  (Fig.  1  and  SI;  Sparks  et  al.,  2013).  These  consistencies 
support  the  assumption  that  putative  chemoreception  gene  tran¬ 
scripts  yielding  positive  RPKM  values  above  1.0  in  labella  and  tarsi 
are  expressed  and  functional  in  these  organs. 

4.  Discussion 

Mosquito  behaviors  such  as  host-seeking,  feeding,  oviposition 
and  mating  are  influenced  by  sensory  inputs  from  labella  and  tarsi 
(Clements,  1992).  These  appendages  harbor  uniporous  cherno- 
sensory  sensilla  as  well  as  sensors  of  physical  stimuli  like  heat  or 
mechanical  force  (Mclverand  Siemicki,  1978;  Lee  and  Craig,  2009). 
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Fig.  3.  Obp  expression  in  labella  and  tarsi.  Numbered  values  represent  the  average  RPKM  values  for  male  and  female  labella  or  male  and  female  tarsi.  Averages  for  tarsi  include  all 
three  tarsal  types:  pro-,  meso-  and  meta-thoracic.  Heat  map  coloration  is  capped  at  30,000  RPKM  with  W%  intensity  at  50  RPKM.  AaegObps  are  divided  into  three  columns 
representing  classic  OBPs,  plus-C  OBPs,  and  two-domain  OBPs.  (*)  indicates  homology  to  D.  melanogaster  LUSH;  (**)  indicates  homology  to  D.  melanogaster  OSE/OSF  (Zhou  et  al, 
2008;  Manoharan  et  al.,  2013). 


4.2.  Olfactory  genes  in  gustatory  organs 

It  remains  unclear  with  which  sensory  neurons  the  AaegORs  and 
AaegIRs  associate  in  the  labellum  and  tarsi  of  mosquitoes.  A  pre¬ 
vious  report  showed  ORCO  immunoreactivity  of  a  few  chemore- 
ceptor  sensilla  on  the  labella  of  Ae.  aegypti  (Melo  et  al.,  2004) 
suggesting  that  AaegORs  function  in  the  context  of  gustation. 
Additionally,  multiple  non-ORCO  Or  transcripts  have  been  ampli¬ 
fied  from  single  gustatory  sensilla  on  the  labella  from  An.  gambiae 
(Kwon  et  al.,  2006),  from  which  at  least  a  few  sensory  neurons 
projected  to  the  antennal  lobe  (Kwon  et  al.,  2006).  Sensilla  on  the 
labellum  of  Ae.  aegypti  and  other  mosquitoes  have  been  described 
as  either  uniporous  or  aporous  (Mclver  and  Siemicki,  1978;  Lee  and 
Craig,  2009),  which  suggests  that  these  volatile-sensing  ORs  and  IRs 
may  function  in  neurons  housed  in  sensilla  that  are  not  multi- 
porous,  as  occurs  on  the  antennae  and  maxillary  palps.  The  sensi¬ 
tivity  of  specific  AaegORs  and  AaegIRs  in  gustatory  hairs  of  the 
labella  and  tarsi  to  volatile  odorants  requires  future  investigation  in 
Ae.  aegypti. 

The  elongated  proboscis  of  mosquitoes  allows  them  to  feed  on 
floral  sugar  sources  as  well  as  vertebrate  hosts  (Clements,  1992). 


Has  expression  of  olfactory  receptors  like  ORs  and  IRs  in  the  labella 
co-evolved  with  the  elongation  of  the  mosquito  proboscis?  It  has 
been  shown  in  An.  maculipennis  that  labella  are  dispensable  for 
normal  blood-feeding  behavior  (Robinson,  1939),  but  the  opposite 
conclusion  was  made  for  Ae.  aegypti  in  a  later  study  (Jones  and  Pilitt, 
1973).  Targeted  knockdown  of  olfactory  receptor  genes  specifically 
in  proboscis  tissue  may  help  elucidate  the  effect  on  behavior  of  ORs 
and  IRs  in  the  labella  of  Ae.  aegypti. 

DEET  has  been  shown  to  block  Ae.  aegypti  host-seeking  behavior 
through  an  ORCO-mediated  pathway  (DeGennaro  et  al.,  2013),  and 
we  have  shown  labellar  expression  of  both  Oreo  and  odorant  acti¬ 
vated  Ors;  therefore,  it  is  possible  there  exists  a  labellar  component 
to  this  ORCO-mediated  avoidance  of  volatile  DEET.  Whether  or  not 
ORs  are  involved  in  the  physiological  response  to  DEET  by  GRNs  of 
the  labella  (Sanford  et  al.,  2013)  requires  direct  investigation.  In 
addition  to  serving  as  an  ion  channel  (Sato  et  al.,  2008;  Wicher 
et  al.,  2008),  ORCO  trafficking  of  conventional  ORs  is  necessary  for 
proper  dendritic  localization  of  ORs  (Larsson  et  al.,  2004).  Perhaps 
ORCO  may  participate  in  the  trafficking  of  other  membrane  pro¬ 
teins  or  function  downstream  of  non-OR  receptor  events  in  the 
labella. 
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Fig.  4.  Chemoreception  gene  enrichment  in  labella  and  tarsi.  RNA  samples  used  for  RT- 
qPCR  were  obtained  from  the  tissues  depicted  in  color.  Each  gene  is  reported  in  'times 
enrichment'  as  a  ratio  of  relative  abundance  of  female  labella  to  female  carcass.  Each 
gene  was  first  normalized  to  housekeeping  gene  AaegLAsP  within  the  given  tissue.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to 
the  web  version  of  this  article.) 


Of  those  species  surveyed  by  intensive  sequence  search  analysis, 
Ae.  aegypti  presents  the  highest  number  of  putative  Irs  (95  genes). 
Only  a  fourth  of  the  Dmellrs  identified  by  sequence  homology  are 
expressed  in  the  antennae  (Rytz  et  al.,  2013),  and  evolutionary 
analysis  of  IRs  throughout  protostomes  suggests  that  the  function 
of  species-specific  IRs  (those  representing  species-specific  gene 
expansions)  extend  beyond  olfaction  to  other  sensory  modalities 
(Croset  et  al.,  2010).  While  IRs  confer  odor  sensitivity  to  specific 
D.  melanogaster  ORNs  and  are  the  defining  molecular  feature  of  a 
large  subset  of  antennal  ORNs  (Benton  et  al.,  2009),  IRs  may  also  be 
sensitive  to  non-volatile  stimuli.  Expression  of  this  ancient  family  of 
receptors  in  chemoreceptive  tissue  extends  to  aquatic  protostomes 
(Croset  et  al.,  2010);  therefore,  IR  involvement  in  the  reception  of 
non-volatile  compounds  very  likely  represents  an  ancestral  func¬ 
tion  of  these  receptors.  Deciphering  whether  the  labellar  or  tarsal 
Aaeglrs  identified  in  our  survey  respond  to  volatile  and/or  non¬ 
volatile  compounds  is  worthy  of  future  consideration.  The  ques¬ 
tion  remains  too  as  to  whether  all  insect  IRs  are  chemoreceptors. 
Thirty  Aaeglrs  show  expression  levels  above  3.0  RPKM  in  male  or 
female  labellar  tissue  samples,  slightly  more  than  the  twenty-four 
AaegGrs  expressing  above  this  threshold  in  our  previous  report 
(Sparks  et  al.,  2013).  Sixteen  Aaeglrs  show  expression  levels  above 
3.0  RPKM  in  at  least  one  male  or  female  tarsal  tissue  sample,  as 
opposed  to  only  eight  AaegGrs  expressing  above  this  threshold 
(Sparks  et  al.,  2013).  Thus,  Irs  may  mediate  many  of  the  gustatory 
signals  affecting  behavior  of  Ae.  aegypti. 

Although  an  IR-mediated  DEET  avoidance  pathway  has  recently 
been  identified  in  the  antennae  of  D.  melanogaster  (Kain  et  al., 
2013),  we  did  not  observe  the  expression  of  the  homolog  Aae- 
glr40a  in  labella  or  tarsi  (Figure  SI).  However,  as  at  least  thirty-six 
putative  Irs  are  expressed  in  the  labella  and/or  tarsi  of  female  Ae. 
aegypti  (Figure  SI),  they  may  be  involved  in  detection  of  repellents 
like  DEET  and  could  provide  potential  targets  for  behavioral 
disruption. 

The  recent  finding  that  OBPs  affect  responses  of  insect  GRNs 
(Jeong  et  al.,  2013;  Swarup  et  al„  2013)  is  not  surprising  given  their 
promiscuous  expression  in  non-olfactory  tissues  (Figs.  1  and  2)  as 
demonstrated  in  our  survey.  Based  on  their  lack  of  expression,  Two- 


Domain  OBPs  do  not  appear  to  affect  gustatory  inputs  in  the  labella 
or  tarsi  adult  Ae.  aegypti  (Fig.  3).  The  Two-Domain  OBP  AaegObp45 
is  expressed  and  accumulates  in  ovaries,  and  is  thought  to  mediate 
eggshell  formation  (Costa-da-Silva  et  al.,  2013).  As  demonstrations 
of  insect  OBP  function  including  specific  ligands  or  other  functional 
data  is  limited  (Grosse-Wilde  et  al„  2006;  Syed  et  al.,  2006;  Forstner 
et  al.,  2009),  AaegOBP  homology  does  not  inform  possible  func¬ 
tionality.  A  few  AaegOBPs  are  homologous  to  well-studied  Dme- 
lOBPs  (LUSH  and  AaegOBP34,  OSE/OSF  and  AaegOBP38,  Manoharan 
et  al.,  2013;  Fig.  3),  but  speculation  as  to  their  function  in  mosquito 
species  requires  further  investigation. 

4.2.  Conclusions 

We  have  revealed  the  expression  of  Ors,  Irs,  Snmp2,  Trp  channels, 
Ppk  channels,  Obps,  CheAs  and  CheBs  in  the  main  gustatory  ap¬ 
pendages  of  Ae.  aegypti,  genes  that  serve  as  useful  targets  towards 
the  behavioral  disruption  of  this  important  disease  vector.  Many 
genes  not  discussed  here  show  low,  but  detectable  levels  of 
expression  (<3.0  RPKM,  Figure  SI).  These  genes  may  be  functional 
at  low  copy  level  or  may  be  active  in  a  relatively  small  subset  of  cells 
present  in  the  tissue  we  surveyed.  Some  chemoreception  genes 
may  be  more  actively  transcribed  at  different  animal  ages  or  life 
stages.  Of  the  genes  assessed  by  RT-qPCR,  all  with  RPKM  values 
greater  than  1  are  highly-enriched  in  labellar  tissue  with  respect  to 
carcass  tissue,  suggesting  that  most  of,  if  not  all,  the  genes  identi¬ 
fied  in  our  survey  are  involved  in  chemoreception.  The  next  steps  in 
characterizing  the  chemoreception  capability  of  the  labella  and 
tarsi  of  Ae.  aegypti  include  determining  the  sensitivities  of  specific 
receptors  (GRs,  ORs,  and  IRs),  assessing  the  relationship  of  these 
receptors  to  putative  ion  channels  like  Trps  and  Ppks,  and  eluci¬ 
dating  the  function  of  highly  expressed  accessory  genes  like  Snmp2 
and  the  most  abundant  Obps  and  CheA/Bs. 
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